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Abstract 

A  performance  model  for  a  reformer,  consisting  of  the  catalytic  partial  oxidation  (CPO),  high-  and  low-temperature  water-gas  shift  (HTWGS 
and  LTWGS),  and  preferential  oxidation  (PROX)  reactors,  has  been  formulated.  The  model  predicts  the  composition  and  temperature  of  the 
hydrogen-rich  reformed  fuel-gas  mixture  needed  for  the  fuel  cell  applications. 

The  mathematical  model  equations,  based  on  the  principles  of  classical  thermodynamics  and  chemical  kinetics,  were  implemented  into 
a  computer  program.  The  resulting  software  was  employed  to  calculate  the  chemical  species  molar  flow  rates  and  the  gas  mixture  stream 
temperature  for  the  steady-state  operation  of  the  reformer.  Typical  computed  results,  such  as  the  gas  mixture  temperature  at  the  CPO  reactor 
exit  and  the  profiles  of  the  fractional  conversion  of  carbon  monoxide,  temperature,  and  mole  fractions  of  the  chemical  species  as  a  function 
of  the  catalyst  weight  in  the  HTWGS,  LTWGS,  and  PROX  reactors,  are  here  presented  at  the  carbon-to-oxygen  atom  ratio  (C/O)  of  1  for  the 
feed  mixture  of  n-decane  (fuel)  and  dry  air  (oxidant). 

©  2004  Elsevier  B  .  V.  All  rights  reserved. 
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1.  Introduction 

Due  to  the  global  enhanced  concern  about  conserving  the 
earth  energy  resources  and  reducing  pollution,  the  fuel  cell  is 
a  promising  technology  to  replace  the  conventional  internal 
combustion  engine  for  power  production.  The  ideal  fuel  for 
efficient  performance  is  pure  hydrogen  in  low  temperature 
fuel  cells  since  it  simplifies  system  integration,  maximizes 
system  efficiency,  and  produces  zero  harmful  emissions  [1]. 
However,  hydrogen  is  not  naturally  available  as  a  fuel,  but  it 
is  stored  in  high  density  in  hydrocarbon  fuels.  Processing  of 
hydrocarbon  fuels  is,  therefore,  necessary  to  extract  hydrogen 
for  stationary  and  mobile  fuel  cell  applications. 

Natural  gas  and  petroleum  liquids  contain  sulfur  com¬ 
pounds  that  are  harmful  to  a  catalyst  used  in  a  fuel  reformer 
unit  to  produce  hydrogen.  Therefore,  sulfur  must  be  removed 
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from  a  sulfur-containing  fuel  to  a  level  tolerable  to  a  fuel  re¬ 
former  catalyst.  Of  the  three  common  methods,  (a)  steam 
reforming  (SR),  (b)  catalytic  partial  oxidation  (CPO),  (c) 
autothermal  or  oxidative  steam  reforming  (ATR)  to  gener¬ 
ate  hydrogen-containing  gas  mixture  for  fuel  cell  applica¬ 
tions,  the  CPO  process  is  less  energy  consuming  compared 
to  other  processes  and  has  simpler  system  configuration.  Cat¬ 
alytic  partial  oxidation  reactions  are  much  faster  than  steam 
reforming  reactions.  Therefore,  smaller  reactors  can  be  re¬ 
alized  with  high  throughput  for  hydrogen  production  [2]. 
Catalytic  partial  oxidation  of  hydrocarbon  fuels  (e.g.  diesel, 
gasoline,  ft-decane,  ^2-hexadecane)  has  been  studied  recently 
[3-5].  The  results  from  the  5wt.%  Rh-on-7-Al203  catalyst 
in  a  short  contact  time  (millisecond  scale)  monolith  reactor, 
have  shown  reaction  pathways  changing  from  combustion 
at  low  carbon  to  oxygen  atom  ratio  (C/O  1)  to  syngas 
production  around  (C/O  ^1),  and  then  to  olefin  production 
at  higher  (C/O)  ratios  [3].  Almost  complete  conversion  of  a 
fuel  and  oxygen  are  achievable  with  selectivity  higher  than 
80%  for  syngas  production  for  (C/O  atom  ratio)  ~  1 .  In  the 
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Nomenclature 

activity  of  chemical  species  i  in  the  gas-solid 
carbon  system  at  the  equilibrium  state 
Ap(i) ,  Z?F(i)  first  and  second  heat  capacity  coefficients  of 
the  fuel  in  the  liquid  phase 

Ai,  Bj ,  Cu  Di  constant  pressure  heat  capacity  coeffi¬ 
cients  of  the  chemical  species  in  the  gas  phase 

ACPO,exit’  fiCK). exit-  ccpo,exit’  DCPO,exit  parameters  as 
defined  in  Section  2,  following  Eq.  (4) 

C°n  •  constant  pressure  standard- state  heat  capacity 

of  a  gaseous  chemical  species  i  at  the  system 

temperature  (Jmol-1  K-1) 

<iint,wGSR  internal  diameter  of  a  water-gas  shift  reactor 

(m) 

Dp  average  diameter  of  particles  in  a  packed-bed 
catalytic  reactor  (m) 

/i,  /2  functions  defined  in  Eqs.  (3c)  and  (3d) 

G  superficial  mass  flow  velocity  in  a  packed-bed 
reactor  (kg  m-2  s-1) 

G  -  standard- state  Gibbs  free  energy  of  a  chemi¬ 
cal  species  i  at  the  system  temperature  (J  (g- 
mol)-1) 

AGg  standard- state  Gibbs  free  energy  change  of  the 
Boudouard  reaction  (J  (g-mol)-1) 

A^wgsr  standard- state  Gibbs  free  energy  change  for 
the  water-gas  shift  reaction  (7)  (J  (g-mol)-1) 
H°T  ,  HfT  standard- state  molar  enthalpy  of  a  gaseous 
chemical  species  i  at  a  reference  temperature, 
To  (e.g.  To  =  298.15  K),  and  at  temperature,  T , 
respectively  (J  (g-mol)-1) 

A  H°xn  standard- state  enthalpy  of  a  reaction  (J(g- 
mol)-1) 

A  //r°xn  CPO  standard- state  enthalpy  of  the  partial  ox¬ 
idation  reaction  (2)  at  the  reference  tempera¬ 
ture,  To  (=  298.15  K)  (J  (g-mol)-1) 

A//°n  wgsr  standard- state  enthalpy  of  the  water-gas 

shift  reaction  (7)  (J  (g-mol)-1) 

A  //r°xn  PROX  standard-state  enthalpy  change  of  the  CO 
oxidation  reaction  (18)  occurring  in  the  prefer¬ 
ential  oxidation  (PROX)  reactor  (J  (g-mol)-1) 
Kq? b  ,  Xe  wGSR  equilibrium  constants  for  the  B  oudouard 
and  water-gas  shift  reactions,  respectively 
Mj  molecular  weight  of  a  chemical  species  j 
(kg  (g-mol)-1) 

Nj  molar  flow  rate  of  a  chemical  species  j  in  a  gas 
mixture  stream  (g-mol  s-1) 

Pi  partial  pressure  of  a  chemical  species  i  in  the 
gas-phase  mixture  (bar) 

P  total  pressure  (bar) 

Pinj/vap  total  pressure  in  the  water  injector/vaporizer 
unit  (bar) 

P°  standard- state  pressure  (=  1  bar) 


2htwgsr-ltwgsr  heal  removal  rate  in  the  external  in¬ 
tercooler  (I)  located  between  the  high-  and  low- 
temperature  water-gas  shift  reactors  (J  s-1) 

Q  ltwgsr-prox  heat  removal  rate  in  the  external  in¬ 
tercooler  (II)  located  between  low-temperature 
water-gas  shift  and  preferential  oxidation  reac¬ 
tors  (J  s-1) 

— rco.HTWGSR  rate  of  CO  conversion  in  the  high- 
temperature  water-gas  shift  reactor  (g-mol  (kg- 
cat.)-1  s-1) 

— rco.LTWGSR  rate  of  CO  conversion  in  the  low- 
temperature  water-gas  shift  reactor  (g-mol  (kg- 
cat.)-1  s-1) 

— rco.TOF  turnover  frequency  for  CO  oxidation  in  the 
preferential  oxidation  reactor  ((g-mol  of  CO) 
(g-mol  of  Pt  atoms)-1  s-1) 

R  gas  constant  (8.314  J  g-mol-1  K-1) 

T,  To  temperature,  the  reference  state  temperature 
(e.g.,  To  =  298.15  K),  respectively  (K) 

Tnbpt.F  normal  boiling  point  of  a  liquid  fuel  (F)  (K) 

U  overall  heat  transfer  coefficient 
(Jm-2  K-1  s-1) 

VW(l)  molar  volume  of  the  liquid  water  (m3  g-mol- 1 ) 
W  weight  of  a  catalyst  in  a  catalytic  reactor  (kg) 
Xco,  Xco,eq  actual  and  equilibrium  fractional  conver¬ 
sion  of  CO,  respectively 

yt  mole  fraction  of  a  chemical  species  in  the  gas 
phase  mixture 

Greek  letters 

cij  ratio  of  molar  flow  rate  of  a  chemical  species  j 
to  the  molar  flow  rate  of  CO 
p  coefficient  of  expansion  of  liquid  water  =  638.6 

X  10-6  K-1  at  100  °C 

Vi  stoichiometric  number  of  a  chemical  species  in 
a  chemical  reaction 

0b  void  fraction  of  the  bed  in  a  catalytic  reactor 
Pb  bulk  density  of  the  catalyst  in  a  packed-bed 

_ o 

reactor  (kg  mb  ) 

pii  viscosity  of  a  pure  component  i  at  the  gas  mix¬ 
ture  temperature  and  pressure  conditions  (Pa  s 
(=  kgm-1  s-1)) 

Pmix  viscosity  of  a  gas-phase  mixture  (Pa  s) 

< pp  a  dimensionless  quantity  defined  in  Eq.  (15c) 


reformer  modeling  work  presented  in  this  paper,  partial  oxi¬ 
dation  of  ft-decane  on  the  Rh-on-y-Al203  catalyst  is  assumed 
to  take  place.  It  is  mentioned  here  that  at  the  time  of  com¬ 
pletion  of  this  work,  the  chemical  kinetics  rate  law  equation 
for  the  partial  oxidation  of  ft-decane  on  this  catalyst  was  not 
available  for  its  application  in  the  performance  model  pre¬ 
sented  in  the  next  section. 
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The  required  degree  of  processing  of  the  hydrogen-rich 
reformate  from  a  catalytic  partial  oxidation  reactor,  e.g.  Rh- 
on-y-Al203  catalytic  reactor,  is  an  issue  related  to  the  type 
of  a  fuel  cell.  High  temperature  fuel  cells,  such  as  solid  oxide 
and  molten  carbonate  fuel  cells,  can  directly  use  the  refor¬ 
mate  from  a  partial  oxidation  reactor.  The  requirement  on 
the  CO  concentration  in  the  fuel  feed  to  the  platinum  anode 
of  the  phosphoric  acid  fuel  cell  (PAFC)  is  0.25-0.5  mol%  of 
CO.  The  CO  tolerance  for  the  low  temperature  proton  ex¬ 
change  membrane  fuel  cell  (PEMFC)  is  even  lower;  about 
10  ppm.  Presence  of  significant  amount  of  CO  in  the  fuel 
feed  to  the  anode  side  of  the  low  temperature  fuel  cells 
(PAFC  and  PEMFC)  increases  the  anode  overpotential;  con¬ 
sequently,  resulting  in  the  reduction  of  the  cell  operational 
voltage  and  the  cell  electric  power  production.  The  high-  and 
low-temperature  stages  of  the  catalytic  CO  conversion  by  its 
reaction  with  steam  via  shift  reaction  are  adequate  to  decrease 
the  carbon  monoxide  content  in  the  reformate  from  a  catalytic 
partial  oxidation  reactor  to  meet  the  requirement  of  the  PAFC. 
Carbon  monoxide  removal  is  required  after  the  shift  reactors 
for  PEM  fuel  cells  by  a  method,  such  as  the  selective  or  pref¬ 
erential  oxidation  reaction.  A  small  amount  of  air  (typically 
about  2  mol%)  [6]  is  added  to  the  fuel  stream  exiting  from 
a  low-temperature,  catalytic  shift  reactor;  which  then  passes 
over  a  precious  metal  catalyst,  such  as  the  Pt-on-Al203  cata¬ 
lyst.  This  catalyst  preferentially  adsorbs  CO,  rather  than  H2. 
Carbon  monoxide,  in  the  chemisorbed  state  on  the  catalytic 
active  sites,  reacts  with  the  oxygen  of  the  added  air  for  its 
conversion  to  CO2. 

The  catalytic  carbon  monoxide  conversion,  via  the  water- 
gas  shift  reaction  (CO  (g)  +  H20  (g)  H2  (g)  +  C02  (g)), 
is  believed  to  take  place  by  only  two  possible  mechanisms. 
In  the  associative  mechanism,  the  reactants  are  chemisorbed 
on  the  active  sites  followed  by  decomposition  to  the  product 
species.  According  to  the  second,  regenerative  mechanism, 
the  cyclic  reduction-oxidation  of  the  catalytic  material  leads 
to  the  formation  of  the  product  species.  Both  mechanisms  are 
supported  for  the  LTWGS  catalysts;  whereas  the  regenerative 
mechanism  is  supported  for  the  HTWGS  catalysts  [7,8]. 

Steam  to  carbon  monoxide  molar  ratio  has  a  critical  impact 
on  the  water-gas  shift  reaction.  For  the  catalytic  HTWGSR 
(e.g.  Fe3C>4  on  Cr203)  operation  under  low  H2O/CO  molar 
ratio,  formation  of  metallic  iron  is  possible.  Consequently, 
methanation,  carbon  deposition,  and  Fischer-Tropsch  pro¬ 
cess  can  occur  [7].  Copper-based  catalysts  (e.g.  Cu-on- 
Zn0/Al203)  for  the  LTWGSR  have  no  tendency  to  catalyze 
the  Fisher-Tropsch  reactions  and  have  very  low  activity  for 
breaking  C— O  or  forming  C— C  bond.  Therefore,  coke  forma¬ 
tion  or  carbon  deposition  on  Cu-based  LTWGSR  catalysts  is 
not  a  serious  problem  [9].  Use  of  a  membrane,  with  simul¬ 
taneous  hydrogen  separation  coupled  with  WGSR,  results  in 
the  enhancement  of  CO  conversion  when  the  H20-to-CO  mo¬ 
lar  ratio  is  low.  In  comparison  with  the  industrial  HTWGS  and 
LTWGS  catalysts,  the  experimental  results  have  shown  that 
Pt-on-Zr02  has  a  potential  for  WGSR  under  low  H20-to-CO 
molar  ratio  [10].  Research  has  been  in  progress  to  explore 


other  effective  catalysts  for  WGSR.  Non-precious  metals, 
such  as  Ni  and  Cu  on  cerium  oxide,  have  shown  higher  activ¬ 
ity  than  the  conventional  LTWGSR  catalyst  (Cu-on-Al203), 
for  example,  at  high  temperatures  [11,12];  whereas  precious 
metals,  such  as  platinum  group  metals  on  cerium  oxides,  have 
shown  activity  in  a  wide  temperature  range  [13].  Nano-scale 
crystals  of  cerium  oxides  with  precious  and  non-precious 
metals  enhance  the  WGSR  rate.  The  promotion  effect  is  at¬ 
tributed  to  the  oxygen  storage  capacity  of  the  oxides,  when 
they  reversibly  exchange  oxygen  ions  in  the  CO  oxidation 
reaction  [14,15].  When  noble  metals  are  used  with  cerium 
oxide,  the  WGS  reactors  can  be  operated  at  higher  temper¬ 
atures  that  are  favorable  to  the  kinetics.  Thus,  the  disadvan¬ 
tages  of  the  conventional  WGS  catalysts  can  be  eliminated. 
However,  currently,  there  are  no  kinetic  rate  law  expressions, 
describing  the  WGSR  on  these  catalysts,  available. 

There  have  been  extensive  efforts  to  find  effective  and 
selective  catalysts  for  preferential  oxidation  of  CO  in  the 
hydrogen-rich  gas  mixtures  [16,17].  The  relevant  informa¬ 
tion  on  the  catalytic  preferential  oxidation  of  CO  over  the 
hydrogen  oxidation  reaction  is  found  in  [18,19,21].  It  has 
been  found  that  the  CO  oxidation  on  Pt-on-Al203  catalyst 
increases  with  an  increase  in  the  reaction  temperature  to  a 
maximum  value,  then  decreases  with  a  similar  trend  of  the 
CO  oxidation  selectivity.  This  behavior  is  explained  as  fol¬ 
lows.  Carbon  monoxide  inhibits  adsorption  of  H2  at  temper¬ 
atures  below  200  °C.  At  temperatures  higher  than  200  °C, 
the  CO  adsorption  is  retarded,  whereas  H2  adsorption  is  en¬ 
hanced.  This  results  in  lowering  of  the  catalyst  CO  oxidation 
selectivity  relative  to  hydrogen  oxidation. 

It  has  been  shown  that  CO  oxidation  is  a  function  of  tem¬ 
perature,  gas  mixture  composition,  total  gas  mixture  flow 
rate,  and  catalyst  loading;  whereas  the  CO  oxidation  selec¬ 
tivity  is  dependent  only  on  the  reactant  feed  gas  composition 
and  temperature.  Methanation  of  CO  and  CO2,  water-gas 
and  reverse  water-gas  shift  reactions  can  possibly  proceed 
on  the  Pt-on-Al203  catalyst.  However,  these  reactions  occur 
insignificantly  in  the  typical  CO  selective  temperature  range 
150-250  °C  [20,21];  a  temperature  range  of  interest  for  the 
modeling  work  reported  in  this  paper. 

The  fuel  reformer  performance  model  reported  here  is  a 
preliminary  attempt  to  quantitatively  simulate  a  fuel  proces¬ 
sor  to  supply  hydrogen-rich  fuel  gas  mixture  to  a  fuel  cell 
stack  operating  at  a  relatively  low  temperature  (e.g.  solid 
polymer  electrolyte  phosphoric  acid  fuel  cell  stack).  Fig.  1 
shows  the  flow  diagram  of  the  process  units  of  the  fuel  re¬ 
former.  A  liquid  fuel  (for  example,  ft-decane)  is  vaporized 
and  heated/mixed  with  air  in  a  heater- vaporizer  unit.  The  gas 
mixture  leaving  the  heater- vaporizer  unit  is  fed  to  the  mono¬ 
lith,  Rh-on-y-Al203,  catalytic  partial  oxidation  microreac¬ 
tor.  Partial  oxidation  of  the  fuel,  ft-decane,  on  the  Rh-on-y- 
AI2O3  catalyst  is  so  fast  that  the  conversion  of  ft-decane  is 
100%  [4]  at  the  carbon-to-oxygen  atom  ratio  of  1.  To  avoid 
sintering  of  the  HTWGSR  catalyst,  Fe304-on-Cr203,  tem¬ 
perature  of  gas  exiting  from  the  CPO  reactor  is  decreased  to 
an  acceptable  level  (suitable  HTWGSR  operational  temper- 
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ature  range:  320-400  °C  [22]).  This  requirement  is  achieved 
by  water  injection,  directly,  into  the  hot  gas  stream  exiting 
from  the  CPO  reactor.  Water  is  also  needed  for  the  CO  con¬ 
version  via  water-gas  shift  reaction.  The  needed  water  may 
be  supplied  from  the  cathode  side  of  a  fuel  cell  stack  being 
fueled  by  the  reformer.  An  external  inter-cooler  (EXT-I)  is 
necessary  between  the  HTWGSR  and  LTWGSR  to  prevent 
sintering  of  the  LTWGSR  catalyst,  Cu-on-ZnO/Al203  (the 
operational  temperature  range:  120-250  °C)  [23].  The  gas 
mixture  stream  leaving  the  LTWGSR  is  cooled  in  the  second 
external  intercooler  (EXT-II)  before  it  enters  the  preferen¬ 
tial  oxidation  (PROX)  reactor.  A  small  amount  of  air  (about 
2  mol%)  is  also  added  to  the  gas  mixture  stream,  being  cooled 
in  the  inter-cooler  (EXT-II),  for  the  preferential  oxidation  of 


2.  Preliminary  reformer  performance  mathematical 
model 

The  first  unit  shown  in  the  sketch  of  the  fuel  cell  reformer 
in  Fig.  1  is  the  fuel  vaporizer-mixer-heater.  A  liquid  hydro¬ 
carbon  fuel,  ft-decane,  is  injected,  vaporized  and  mixed  with 
dry  air.  The  fuel/air  mixer  is  heated  to  a  required  inlet  temper¬ 
ature  of  the  Rh-on-alumina  catalytic  partial  oxidation  (CPO) 
reactor.  By  the  use  of  the  first-law  of  thermodynamics,  with 
the  assumptions  of  negligible  potential  and  kinetic  energy 
changes  and  ideal  gas  behavior,  and  no  shaft  work  exchange 
between  the  system  (i.e.  vaporizer-mixer-heater  unit) 
and  surroundings,  the  following  equation  for  the  required 
rate  of  heat  supply  to  the  vaporizer-mixer-heater  unit  was 
derived: 


•  • 
2heat/vap  =  ^F,in 


R 


^F(g)(7cPO,in  -  ?nbpt,F)  + 


#F(g)(^CPO,in  ^nbpt,F^  ,  ^F(g)^CPO,in  ^nbpbF^ 


+ 


+ 


^F(g)(^cPO,in  ^nbpt,F^ 


A  K 


o 


vap,F,rnbpt 


+ 


^F(l)(7nbpt,F  -  ^F(l),in)  + 


#F(l)(^nbpt,F  ^F(l),in) 


+  Miir,CPO,in 


R 


^air(7cPO  ,in  ^air  ,in  )  + 


#air(^CPO,in  ^air,in^  ,  ^air(^CPO,in  ^air.in) 


—  D 


1 


1 


air 


7cPO,in  Tc 


air,m 


+ 


CO  on  the  Pt-on-Al203  catalyst;  with  the  operational  tem¬ 
perature  range:  150-350  °C  [21].  The  product  gas  mixture 
from  the  preferential  oxidation  reactor  may  be  fed  directly  to 
a  fuel  cell  stack,  further  treated  for  obtaining  pure  hydrogen 
(if  needed)  or  stored  in  cylinders  after  cooling  for  its  later  use. 


The  first  term  on  the  right-hand  side  of  Eq.  (1)  repre¬ 
sents  the  amount  of  heat  needed  to  vaporize  and  heat  a  liquid 
fuel,  such  as  ft-decane,  to  a  required  inlet  temperature  of  the 
CPO  reactor.  The  second  term  on  the  right-hand  side  repre¬ 
sents  the  amount  of  heat  needed  to  heat  the  entering  air  to 
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Fig.  1.  Schematic  diagram  of  a  reformer  to  produce  fuel  for  a  fuel  cell. 
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the  CPO  reactor  required  inlet  temperature.  The  fuel  used 
in  our  modeling/numerical  simulation  is  ^-decane.  ft-Decane 
in  gas  phase  reacts  with  oxygen  of  the  oxidant,  air,  accord¬ 
ing  to  the  following  partial  oxidation  reaction  stoichiometric 
equation: 

C10H22  (g)  +  502  (g)  ->  IOCO  (g)  +  1 1H2  (g)  (2) 

Experimental  work  [3]  has  shown  that  the  fuel  and  oxygen 
are  almost  completely  consumed  in  the  Rh-on-Al203  cat¬ 
alytic  partial  oxidation  reactor  at  the  carbon-to-oxygen  atom 


^CPO  ,exit  —  1  + 

where 


A^rxn,CPO,T0  +  (rCPO,in  l)/l 

Jl 


^CPO,in 


^CPO  ,exit 


/  ?CPO, 


exit 


V  T0 


(3a) 
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+  7o(TCPO,in  +  TCPO,in  +  1)  ' 


Cf  +  (5Co?  +  18.81Cn,) 


’  +  (TCPO,in  +  l)(TCPO,in  +  1) 


x  ,  ^^f+5Do2  +  18.81DN2 
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(3c) 


h  =  R 
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A)^CPO,exit 


(3d) 


ratio,  C/O,  of  1.  During  the  period  of  the  modeling  study 
reported  here  the  kinetic  rate  law  equation  for  the  partial  ox¬ 
idation  reaction  (2)  was,  unfortunately,  not  available  [3,25]. 
Therefore,  the  assumption  of  complete  conversion  of  fuel  ( n - 
decane)  to  syngas  production  was  employed  in  the  modeling 
work  presented  here.  By  the  application  of  the  first-law  of 
thermodynamics,  with  the  assumptions  of  negligible  kinetic 
and  potential  energy  changes,  and  idea-gas  behavior;  no  shaft 
work  exchange  between  the  CPO  reactor  system  and  its  sur¬ 
roundings,  and  under  adiabatic  and  steady- state  conditions, 
the  following  mathematical  equation  (3)  was  obtained  to  cal¬ 
culate  the  exit  temperature  of  the  gas  mixture  leaving  the 
CPO  reactor.  In  the  formulation  of  Eq.  (3),  the  air  feed  was 
assumed  to  be  at  the  stoichiometric  level  corresponding  to 
the  carbon-to-oxygen  atom  ratio,  C/O,  of  1  in  the  reaction 
stoichiometric  equation  (2): 


The  Newton-Raphson  method  was  used  to  solve  the  nonlin¬ 
ear  algebraic  equation  (3)  for  Tcpo,exit>  and  hence,  the  hot  gas 
stream  temperature  at  the  CPO  reactor  exit,  7cpo,exit- 

The  hot  gas  mixture  leaving  the  CPO  reactor  is  cooled 
in  the  water  injector/cooler  unit  shown  in  Fig.  1,  by  the  di¬ 
rect  injection  of  liquid  water  into  the  hot  gas  stream  for  ef¬ 
ficient  vaporization  of  liquid  water  and  cooling  of  the  hot 
gas  stream  to  a  desired  inlet  temperature  of  the  HTWGSR 
for  the  catalyst  sintering-free  operation  of  the  reactor.  The 
injected  water  in  the  liquid  phase,  not  only,  meets  the  gas 
stream  cooling  requirement;  but  also,  is  used  to  satisfy  the 
amount  of  water  needed  to  convert  CO  to  CO2  in  the  highl¬ 
and  low-temperature  water-gas  shift  reactors.  Eq.  (4)  was  de¬ 
rived  to  calculate  the  required  amount  of  liquid  water  by  the 
application  of  first-law  of  thermodynamics,  with  the  assump¬ 
tions  of  the  idea  gas  mixture  behavior  and  negligible  kinetic 
and  potential  energy  changes,  with  no  shaft  work  exchange 
between  the  system  (water  injector/cooler)  and  its  surround¬ 
ings  at  the  steady-state,  adiabatic  operating  conditions: 
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where 


is  given  by 


«W,inj 
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The  model  equations  for  the  high-  and  low-temperature 
water-gas  shift  reactors  are  identical  except  for  the  carbon 
monoxide  rate  law  equation  for  its  consumption.  It  is  here 
assumed  that  the  HTWGS  and  LTWGS  packed-bed  catalytic 
reactors  have  plug  flow.  The  mole  balance  on  CO,  over  a  spa¬ 
tial  element  in  the  packed-bed  reactor,  was  applied,  under  the 
steady- state  condition,  to  obtain 


A(7wgsr  —  E  ViG°  O'  =  CO,  H2,  H20,  C02)  (7b) 

i 

where  vCo  =  -1,  vco2  =  1,  vn2o  =  -1,  m2  =  1- 

The  rate  of  CO  disappearance  per  unit  mass  of  the 
Cu/Zn0/Al203  (Cu/Zn  g-atom  ratio  =  2)  catalyst  [23]  in  the 
LTWGSR  is  given  by 


-TCO, LTWGSR 


=  8.22  x  104  exp 


-5701.23 


T 


PCOPH2o(l  -  0 


Thermal  energy  balance  applied  over  a  spatial  element  in  a 
catalytic  bed  water-gas  shift  reactor,  under  the  steady- state 
conditions,  leads  to 


dT 

dW 
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where 

AHvxn  =  J2ViHi  0=CO,H2,H2O,CO2)  (9a) 

i 

U  is  the  overall  heat  transfer  coefficient  for  heat  exchange 
between  the  WGSR  bed  and  its  surroundings. 

Under  the  adiabatic  condition  (i.e.  U  =  0.0),  Eq.  (9)  be¬ 
comes 
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dW 
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The  rate  of  the  CO  disappearance  per  unit  mass  of  the 
Fe304-Cr203  (90-95%  Fe304,  5-10%  Cr203)  catalyst  [22] 
in  the  HTWGSR  is  given  by 


~r  CO,HTWGSR 


=  3.428  x  109  exp 


-11426.51 

T 


cCOcH2o(l  O  (6) 


where  the  reversibility  factor,  f ,  for  the  water-gas  shift  reac¬ 
tion  (WGSR)  is  given  by 


^co2  Pn2  /  Pco  Pr2o 

^e,WGSR 


The  equilibrium  constant,  ^e,WGSR>  for  the  water-gas  shift 
reaction 


CO(g)  +  H20(g)  ^  H2(g)  +  C02(g)  (7) 


Eq.  (10)  predicts  the  gas-phase  temperature  as  a  function  of 
the  catalyst  weight  in  a  water-gas  shift  reactor  bed  under 
the  steady-state,  adiabatic  condition.  The  enthalpy  or  heat 
of  the  water-gas  shift  reaction,  A  H°xn  WGSR,  Is  given  by  the 
following  derived  equation: 

AHr°xn  wgsr  =  -49222.158  +  15.4647  -  0.0022457’2 

9.6775  x  105 

+ - j -  (ID 

It  is  noted  here  that  the  differential  equations  (5)  and  (10), 
should  be  solved  simultaneously  to  calculate  the  CO  actual 
conversion,  Xco,  and  the  bed  gas  temperature,  T,  as  a  function 
of  the  reactor  bed  catalyst  weight.  Focally,  the  gas  phase  and 
the  catalytic  particles  in  the  reactor  bed  are  assumed  to  be  in 
thermal  equilibrium.  Furthermore,  gradients  of  temperature 
and  chemical  species  concentration  in  the  radial  direction  in 
the  tubular  reactor  are  assumed  to  be  zero. 

The  following  expression  for  the  equilibrium  CO  conver¬ 
sion,  Xco,eq>  via  the  water-gas  shift  reaction,  was  derived, 
based  on  the  equilibrium  thermodynamic  fundamentals,  to 
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compare  the  local  actual  and  equilibrium  CO  conversion  val¬ 
ues  in  a  catalytic  water-gas  shift  reactor: 

(1  +  Q'H2o)^e,WGSR  +  (^C02  +  aU2) 


[(1  +  0'H2o)^e,WGSR  +  (o'C02  +  Oiu2)f  ~  4(^e 


AC0,eq  — 

2(/Gj,wgsr  -  1) 
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OiU2  =  . 

AcO,WGSR,in 

Afco,WGSR,in  1 

Oi  co  =  •  =1 

NcO,WGSR,in 

Of  the  two  values  of  Xco,eq  provided  by  Eq.  (12),  only 
the  positive  number  value  <1.0  is  the  acceptable  one.  The 
derived  equilibrium  criterion  for  the  solid  carbon  forma¬ 
tion/deposition  on  the  catalyst  in  the  high-  or  low-temperature 
water-gas  shift  reactor,  via  the  Boudouard  reaction  (13)  is: 

2CO  (g)  ^  C02  (g)  +  C  (s)  (13) 


At  the  gas  mixture  equilibrium  state  with  respect  to  the 
Boudouard  reaction, 


=  exp 


RT  ) 


(i  =  CO  (g),  C02  (g),  C  (s)), 


VCO,B  =  —2,  Vco2,B  =  1,  I'C.B  =  1 


(14) 


AGb  =  X>,bG?  (14a) 

i 

The  activity  of  solid  carbon  is  ac(s )  =  1  and  the  activity  of 
gaseous  species  in  the  ideal  gas-phase  mixture  at  a  pressure, 
such  as  1  bar,  is  given  by 

«(,e  =  ^  0=CO,CO2)  (14b) 

Substitution  for  a^e  from  Eq.  (14b)  into  Eq.  (14),  on  simpli¬ 
fication,  leads  to: 

yco2 

_ylo 

If  the  actual  value  of  [yco2/yco^ actual  §reater  than 
[yco2/yco\q’  there  is  a  possibility  of  a  solid  carbon  forma¬ 
tion/deposition  on  the  catalyst  in  a  water-gas  shift  reactor  ac¬ 
cording  to  the  Boudouard  reaction  [6,26].  If  [yco2/};co] actual 


=  K 
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is  less  than  [ycc^/Acoleq’  no  sohh  carbon  presence  is  pre¬ 
dicted. 


-  t)(^e,WGSRQ'H20  “  Q'C02Q'H2) 

-  (12) 


Total  pressure  drop  can  occur  along  the  length  of  a  water- 
gas  shift  reactor,  predominantly,  due  to  the  presence  of  cat¬ 
alytic  particles  in  the  reactor  bed.  The  well-known  Ergun 
equation  [27]  can  be  used  to  determine  the  total  pressure  as 
a  function  of  the  catalyst  weight,  W,  in  the  reactor: 


d  P 
dW 
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ZX 
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For  the  isothermal  operation  of  a  water-gas  shift  reactor,  with 
negligible  effect  of  pressure  and  composition  on  the  gas  mix¬ 
ture  viscosity;  Eq.  (15)  can  be  integrated  to  obtain 


\  / 1  —  0b\  / 150/4(1  ~0b) 

A  K  A  A> 


-1 0.5 


1.75G  W 


(15a) 


The  gas  mixture  viscosity,  /im ix,  if  not  provided,  can  be  cal¬ 
culated  from  the  following  equation  [28]: 


(15b) 


where 


(15c) 


It  should  be  noted  that  a  pure  component  viscosity,  /x;,  is 
at  the  gas  mixture  temperature  and  pressure  conditions. 
Furthermore,  it  has  been  shown  that  the  semiempirical 
formula,  Eq.  (15b),  reproduces  measured  values  of  the 
viscosities  of  the  gas  mixtures  within  an  average  deviation  of 
about  2%. 

External  intercoolers  (EXT-I  and  EXT- II,  see  Fig.  1)  are 
needed  between  the  HTWGS  and  LTWGS  and  between  the 
LTWGS  and  PROX  reactors  to  adjust  the  reactive  mixture 
stream  temperature  to  the  desired  levels  prior  to  its  entry  to 
the  LTWGS  and  PROX  reactors.  One  should  make  it  cer¬ 
tain  that  the  reactive  gas  mixture  temperatures  at  the  inlets 
of  the  LTWGS  and  PROX  reactors  are  such  that  the  tem¬ 
perature  profiles  in  these  reactors  are  within  the  operational 
temperature  ranges  of  the  catalysts  in  these  reactors  to  prevent 
sintering  of  the  catalyst  and  its  activity  decay.  By  the  applica¬ 
tion  of  first-law  of  thermodynamics;  with  the  assumptions  of 
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negligible  kinetic  and  potential  energy  changes,  and  ideal  gas 
mixture  behavior;  and  no  shaft  work,  the  following  equation 
has  been  derived  for  the  required  heat  removal  rate  in  the  ex¬ 
ternal  intercooler,  EXT-I,  between  the  HTWGS  and  LTWGS 
reactors,  for  the  cooler  steady- state  operation: 

nEXT-l 

^  HTWGSR— LTWGS  R 


The  preferential  oxidation  reactor  is  used  to  decrease  the 
CO  concentration  in  the  hydrogen-containing  gas  stream  to 
an  acceptable  level  for  its  ultimate  application  in  a  low- 
temperature  fuel  cell  stack  (for  example,  phosphoric  acid 
electrolyte  fuel  cell  stack).  Carbon  monoxide  oxidation  takes 
place  on  the  Pt-on-Al203  catalyst  via  the  following  reaction 
stoichiometric  equation: 
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X  ^LTWGSR,  in  ^HTWGSR,  exit  + 1)  + 


J2  jN j,  HTWGSR, exit  C  j 


3  2 

^LTWGSR,  in  (  rHTWGSR,  exit 


+  ^HTWGSR  ,exit  +  D 


+ 


EAiDj 


^LTWGSR,in  ^HTWGSR,  exit 


(16) 


(j  =  H2,  H2O,  CO,  CO2,  N2,  and  O2  if  left  unreacted  in  the 
CPO  reactor), 
where 

THTWGSR,exit  ...  x 

^HTWGSR  ,exit  —  ^  (16a) 

2  LTWGSR,  in 

Eq.  (16)  is  a  simple  algebraic  equation,  the  required  heat 
removal  rate  in  the  external  intercooler,  EXT-I,  can  be  calcu¬ 
lated  by  simply  substituting  the  values  of  the  chemical  species 
molar  flow  rates,  exit  temperature  of  the  reactive  mixture 
stream  leaving  the  HTWGSR,  its  required  temperature  at  the 
inlet  of  the  LTWGSR,  and  the  heat  capacity  coefficients  for 
the  chemical  species  present  in  the  gas  mixture  stream. 

If  the  gas  mixture  stream  leaving  the  LTWGSR  does  not 
contain  sufficient  amount  of  oxygen  to  convert  CO  in  the 
PROX  reactor  to  a  desired  conversion  level,  oxygen  or  air  is 
added,  at  a  controlled  rate,  to  the  hot  gas  stream  either  prior 
to  the  external  intercooler,  EXT-II,  or  in  concurrence  with 
its  cooling  in  the  cooler.  Addition  of  oxidant  in  this  manner 
should  result  in  the  lesser  required  heat  removal  rate  in  the 
EXT-II  cooler.  Air  may  be  added  at  its  ambient  temperature. 
The  derived  equation  for  the  heat  removal  rate  in  the  external 
intercooler,  EXT-II,  for  its  steady-state  operation  is: 


C0(g)  +  ±02(g)^C02(g)  (18) 

This  reaction  (18)  is  selectively  accelerated  relative  to  the 
hydrogen  oxidation  reaction  on  this  catalyst. 

H2(g)+i02(g)4»H20(g)  (18a) 

The  intrinsic  chemical  reaction  rate  law  equation  for  CO  ox¬ 
idation  [21]  is: 

-rco.TOF  =  3.1  X  108  exp  (  9^2'6^  Pco5 Po/  (19> 

This  rate  law  equation  is  for  the  Pt-on-Al203  catalyst,  with 
10%  (by  weight)  of  the  platinum  loading  and  60%  dispersion 
of  platinum  atoms.  The  above  kinetic  rate  law  equation  is 
valid  over  the  temperature  range:  423-623  K.  The  plug  flow 
type,  carbon  monoxide  mole  balance  and  thermal  energy  bal¬ 
ance  equations  similar  to  Eqs.  (5)  and  (10),  respectively,  are 
used  to  predict  the  CO  conversion  and  the  gas  temperature  as  a 
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where 
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Cpi  d  T 


(17a) 


The  first  and  second  sum  terms  within  the  square  brackets 
on  the  right-hand  side  of  Eq.  (17)  are  the  enthalpy  flow  rates 
associated  with  the  hot  gas  stream  leaving  the  LTWGSR  and 
the  oxidant  stream,  respectively.  The  third  sum  term  is  the 
enthalpy  flow  rate  associated  with  the  stream  entering  the 
PROX  reactor. 


function  of  the  catalyst  weight  in  the  PROX  reactor.  Locally, 
thermal  equilibrium  between  the  gas  and  catalytic  particle 
solid  phase  is  assumed.  The  enthalpy  or  heat  of  reaction  (18), 
A//r°xn  PROX,  is  given  below  for  its  use  in  the  thermal  energy 
equation  for  the  PROX  reactor: 


A//r°n  prox  =  -286543.83  +  2.1747  +  9.77  x  10"472 


+ 


8.42  x  105 
T 


(20) 
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3.  Numerical  simulation 

A  computer  program  was  developed  in  Microsoft  Visual 
Basic  6.0  to  generate  numerical  data  from  the  model  perfor¬ 
mance  equations  [29].  The  input  parameters  to  the  computer 
program  are  the  inlet  molar  flow  rates  of  the  fuel,  ft-decane, 
and  air;  the  system  pressure,  inlet  temperature  of  the  fuel  and 
air  to  the  heater-vaporizer  unit,  inlet  temperature  of  the  gas 
mixture  to  the  catalytic  partial  oxidation  (CPO)  reactor,  tem¬ 
peratures  at  the  inlets  of  the  HTWGS,  LTWGS  and  PROX 
reactors,  and  molar  flow  rate  ratio,  O2/CO,  at  the  inlet  of 


the  PROX  reactor.  It  is  here  mentioned  that  some  of  the  re¬ 
former  performance  model  equations  presented  in  Section  2, 
for  example,  the  solid  carbon  presence  criterion  and  reac¬ 
tor  pressure  drop  equations,  were  not  in  the  final  computer 
simulation  stage  at  the  time  of  writing  of  this  paper.  Also, 
the  reference  [24,26,29]  should  be  cited  for  other  parameter 
numerical  data. 

The  heater-vaporizer,  water  injector/cooler,  and  external 
cooler  EXT-I  and  EXT-II  model  equations  are  simple  alge¬ 
braic  equations.  The  solution  of  these  equations  is  straightfor¬ 
ward.  Newton-Raphson  method  is  used  to  estimate  the  exit 


Start 


Inlet  n-decane  fuel  molar  flow  rate,  C/O  atom  ratio  to  (CPO), 
02/C0  molar  ratio  to  (PROX),  inlet  temperatures  for  all  system 

units. 


Set  heat  capacity  coefficients,  molecular  weight, 
standard  heat  of  formation  of  chemical  species,  heat  of 
vaporization  of  the  liquid  n-decane  fuel  and  water. 


Calculate  the  exit  molar  flow  rate  of  every  species  from  CPO  and  setT  cp0 


exit 


=  \ 


Calculate  a  new  value  forT 


CPO  ,  exit  ‘ 


Compare  the 
old  value  with 
the  new  one. 


Set  the  new  value  of 

^ CPO, exit  as  an  old 
(assumed)  value. 


Print  in  the  summary  file  the  CPO  exit 
temperature,  the  exit  molar  flow  rate  of 
every  chemical  species  from  CPO. 

T 

Calculate  the  required  amount  of  liquid  water  needed  by 
the  injector/cooler  unit. 


Print  in  the  summary  file 
the  required  amount  of 
water  needed  for  the 
injector/cooler  unit. 

~  ▼ 

Calculate  the  thermodynamic  CO  equilibrium  conversion,  and 
check  the  carbon  deposition  criterion  at  the  inlet  of  the  HTWGSR. 


Print  in  HTWGSR  file  catalyst  weight,  the  CO 
equilibrium  conversion,  carbon  deposition 
criterion  result,  and  the  mole  fractions  of  all 
chemical  species  at  the  inlet  of  HTWGSR. 


1 


Fig.  2.  Flow  diagram  of  the  program  steps  for  the  fuel  reformer  modeling. 
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Calculate  the  actual  CO  conversion,  temperature,  mole 
fractions  of  chemical  species  and  the  equilibrium  CO 
conversion  at  the  new  catalyst  weight  in  HTWGSR. 


Calculate  the  required  cooling  duty  by  the 
cooler  after  HTWGSR  for  the  gas  mixture. 


Print  in  the  summary  file  the 
required  cooling  duty. 


Calculate  the  equilibrium  CO  conversion 
and  check  the  carbon  deposition  criterion 
at  the  inlet  of  the  LTWGSR. 


Print  in  LTWGSR  file  the 
catalyst  weight,  the  CO 
equilibrium  conversion,  the 
actual  CO  conversion,  carbon 
deposition  criterion,  the  mole 
fraction  of  every  specie. 


Set  the  new  calculated 
actual  CO  conversion  and 
temperature  values  as  old 
values  for  the  new  loop. 


Print  in  HTWGSR  file,  the 
equilibrium  CO  conversion,  the 
actual  CO  conversion,  carbon 
deposition  criterion, 
temperature,  chemical  species 
mole  fractions  at  the  new 


Set  the  new 
calculated  actual 
CO  conversion  and 
temperature  values 
as  old  values  for  the 
new  loop. 


Calculate  the  actual  CO  conversion,  temperature,  equilibrium 
CO  conversion,  and  mole  fractions  of  every  chemical  species  at 
the  new  catalyst  weight  in  LTWGSR. 


Print  in  LTWGSR  file  the 
CO  equilibrium  conversion, 
actual  CO  conversion, 
temperature,  and  the  mole 
fractions  at  the  new  catalyst 
weight. 


Fig.  2.  ( Continued ) 


temperature  from  the  CPO  reactor.  The  coupled  CO  mole  bal¬ 
ance  and  energy  equations  are  solved  simultaneously  by  the 
fourth-order  Runge-Kutta  method  in  the  HTWGS,  LTWGS, 
and  PROX  reactors.  Fig.  2  shows  the  flowchart  of  the  execu¬ 
tion  steps  in  the  program. 


4.  Results  and  discussion 

Fig.  3  shows  the  effect  of  the  CPO  reactor  inlet  tempera¬ 
ture  on  the  gas  mixture  temperature  leaving  the  CPO  reactor 


at  the  carbon  to  oxygen  atom  ratio  of  1 .  The  nature  of  the  gas 
mixture  outlet  temperature  versus  the  gas  mixture  inlet  tem¬ 
perature  plot  is  totally  dependent  on  Eq.  (3)  given  in  Section 
2.  It  is  almost  straight  line  over  the  inlet  temperature  range 
480-520  K.  Results  presented  in  Figs.  4  and  5  correspond  to 
the  parametric  conditions  given  in  Table  1 . 

Fig.  4  shows  the  change  of  actual  CO  conversion,  equilib¬ 
rium  CO  conversion,  and  gas  mixture  temperature  as  a  func¬ 
tion  of  the  catalyst  weight  in  the  HTWGSR.  The  gas  mixture 
temperature  increases  with  an  increase  in  the  catalyst  weight 
due  to  the  exothermic  nature  of  the  water-gas  shift  reaction; 
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Fig.  2.  ( Continued ). 


Fig.  3.  The  effect  of  the  inlet  temperature  of  the  gas  mixture  on  the  exit 
temperature  of  the  CPO  reactor. 

Table  1 

Parametric  data  at  carbon-to-oxygen  atom  ratio  of  1 

Fuel  OZ-C10H22)  feed  rate  (g-mols_1) 

Air  feed  rate  (g-mol  s-1) 

System  pressure  (bar) 

Fuel  and  air  inlet  temperature  (°C) 

Temperatures  at  the  inlets  of  CPO,  HTWGS,  LTWGS 
and  PROX  reactors  (°C) 


whereas  the  equilibrium  CO  conversion  decreases  and  the 
actual  CO  conversion  increases  with  the  increase  in  catalyst 
weight.  Since  the  water-gas  shift  reaction  is  reversible  and  the 
actual  CO  conversion  is  limited  by  the  principles  of  equilib¬ 
rium  thermodynamics;  the  point,  where  the  actual  and  equi¬ 
librium  CO  conversion  curves  meet,  is  the  point  of  maximum 
actual  CO  conversion. 

Eq.  Xco 

800 
600 
400  ^ 

200 
0 

0.0000  0.0002  0.0004  0.0006  0.0008 

W,kg 

Fig.  4.  Temperature,  actual  and  equilibrium  CO  conversion  vs.  catalyst 
weight  in  the  HTWGSR. 
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Fig.  5.  Mole  fraction  profiles  of  chemical  species  in  the  gas  mixture  in  the 
HTWGSR  vs.  the  catalyst  weight. 


Fig.  6.  The  effect  of  total  inlet  molar  flow  rate  on  CO  conversion  in  the 
HTWGSR. 

Fig.  5  shows  the  species  mole  fraction  profiles  as  a  func¬ 
tion  of  the  catalyst  weight  in  the  HTWGSR.  Nitrogen  and 
oxygen  (oxygen  in  trace  amount)  mole  fractions  are  constant 
with  respect  to  the  catalyst  weight  along  the  reactor  length; 
because,  they  are  not  involved  in  the  water-gas  shift  reaction 
and  the  total  gas  mixture  molar  flow  rate  does  not  change  due 
to  the  nature  of  the  reaction  stoichiometric  equation.  Mole 
fractions  of  the  product  species,  hydrogen  and  carbon  diox¬ 
ide,  increase;  whereas  those  of  the  reactant  species,  carbon 
monoxide  and  water,  decrease.  Fig.  6  shows  the  effect  of  the 
total  fuel  and  air  molar  feed  rate  on  the  actual  CO  conver¬ 
sion  versus  the  catalyst  weight  profile  in  the  HTWGSR  at  the 
carbon-to-oxygen  atom  ratio  of  1 .  The  total  molar  feed  rate 
setting  numbers  I,  II,  III  and  IV  correspond  to  the  molar  flow 
rates  given  in  Table  2  and  the  other  parametric  data  are  the 
same  as  shown  in  Table  1 . 

From  Fig.  6,  it  is  apparent  that  a  higher  amount  of  the 
catalyst  is  needed  for  a  higher  total  molar  feed  rate  compared 
to  that  for  a  lower  total  molar  feed  rate  to  achieve  the  same 
fractional  conversion  of  CO. 

Figs.  7-9  show  the  plotted,  computed  data  for  the 
LTWGSR  similar  to  the  HTWGSR  data.  The  explanation 
for  these  plots  is  exactly  similar  to  that  given  above  for 
the  HTWGSR  plots.  However,  the  catalyst  weight  in  the 
LTWGSR  is  in  the  order  of  kg  whereas  that  is  in  the  or- 


Fig.  7.  Temperature,  actual  and  equilibrium  CO  conversion  vs.  catalyst 
weight  in  the  LTWGSR. 


Fig.  8.  Mole  fraction  of  the  chemical  species  vs.  catalyst  weight  in  the 
LTWGSR. 


Fig.  9.  Effect  of  the  total  inlet  molar  feed  rate  on  the  profile  of  the  actual 
CO  conversion  vs.  the  catalyst  weight  in  the  LTWGSR. 

der  of  mg  in  the  HTWGSR.  This  is  due  to  the  lower  chemical 
kinetics  rate  affected  by  the  lower  operational  temperature 
range  that  prevails  in  the  catalytic  LTWGS  reactor  bed.  One 
should  carefully  control  the  temperature  of  the  gas  mixture 
entering  the  LTWGSR  so  that  the  temperature  of  the  gas  mix¬ 
ture  leaving  it  is  less  than  260  °C  to  prevent  sintering  of  the 
copper  catalyst.  Also,  the  copper  catalysts  do  not  have  strong 
activity  toward  breaking  of  C— O  bond;  and  hence,  the  danger 
of  carbon  formation/deposition  via  the  Boudouard  reaction 
is  less  compared  to  the  iron-based  HTWGSR  catalyst  [7]. 

The  preferential  oxidation  (PROX)  reactor  results  shown 
in  Figs.  10-13  correspond  to  the  reformer  parametric  feed 


Table  2 


Reactant  feed  data  at  the  carbon-to-oxygen  atom  (C/O)  ratio  =  1 


Reactant  feed  rate 
setting  number 

Fuel  (n-CioH22)  molar 
flow  rate  (g-mols-1)  (A) 

Air  feed  flow  rate 
(g-mols 1)  (B) 

Total  molar  feed 

rate  =  A  +  B  (g-mols 1) 

I 

2.717E— 04 

6.4696E— 03 

6.7413E— 03 

II 

1.087E— 03 

2.5879E— 02 

2.6966E— 02 

III 

3.261E— 03 

7.7636E— 02 

8.0897E— 02 

IV 

6.0E-03 

14.2857E— 02 

14.8857E— 02 
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Fig.  10.  CO  conversion  and  temperature  profiles  vs.  catalyst  weight  in  the 
PROX  reactor. 


Fig.  11.  Mole  fraction  of  the  chemical  species  vs.  catalyst  weight  in  the 
PROX  reactor. 

rate  data  given  in  Table  1.  Fig.  10  shows  the  actual  CO  con¬ 
version  and  temperature  profiles  as  a  function  of  the  catalyst 
weight  in  the  PROX  reactor  for  the  inlet  oxygen  to  carbon 
monoxide  molar  ratio  of  2.  Close  to  the  reactor  inlet,  the  CO 
conversion  and  the  gas  mixture  temperature  increase  almost 
linearly  as  a  function  of  the  Pt-on-Al203  catalyst  weight.  As 
the  CO  concentration  decreases  with  an  increase  in  the  cat¬ 
alyst  weight,  the  rise  in  the  actual  CO  conversion  becomes 
relatively  steep.  The  gas  mixture  temperature  follows  a  sim¬ 
ilar  trend.  This  observed  behavior  is  in  agreement  with  [30]. 
The  CO  oxidation  rate  has  a  negative  order  with  respect  to 


Fig.  12.  Inlet  gas  mixture  temperature  effect  on  CO  conversion  in  the  PROX 
reactor. 


Fig.  13.  Effect  of  the  inlet  molar  O2/CO  ratio  on  CO  conversion  in  the  PROX 
reactor. 


CO  concentration.  An  increase  in  the  bed  temperature  and  a 
decrease  in  the  CO  concentration  in  the  ppm  range  result  in 
a  relatively  high  value  of  the  CO  oxidation  rate,  and  hence, 
in  the  actual  CO  conversion. 

The  profiles  of  the  mole  fractions  of  chemical  species  as 
a  function  of  the  catalyst  weight  in  the  (PROX)  reactor  are 
shown  in  Fig.  11.  Carbon  monoxide  and  oxygen  being  the 
reactant  species,  their  mole  fractions  decrease  with  an  in¬ 
crease  in  the  catalyst  weight;  whereas  the  mole  fraction  of 
the  product  CO2  increases  with  the  catalyst  weight  increase. 
Hydrogen,  nitrogen,  and  water  species  were  assumed  not  to 
take  part  in  any  reaction.  Therefore,  their  molar  flow  rates  do 
not  change  with  an  increase  in  the  catalyst.  However,  their 
mole  fractions  increase  slightly  with  an  increase  in  the  cata¬ 
lyst  weight  due  to  the  decrease  in  the  total  mixture  molar  flow 
rate  associated  with  the  occurrence  of  the  CO  oxidation  reac¬ 
tion.  Fig.  12  shows  the  effect  of  the  PROX  reactor  inlet  tem¬ 
perature  on  the  actual  CO  conversion  versus  catalyst  weight 
profile.  One  observes  that  the  weight  of  the  catalyst  needed 
is  less  for  a  higher  inlet  temperature  for  the  same  actual  CO 
conversion.  However,  with  an  increase  in  the  bed  tempera¬ 
ture,  CO  desorption  from  the  catalytic  surface  is  enhanced. 
This  affects  the  catalyst  selectivity  towards  the  CO  oxidation 
relative  to  the  hydrogen  oxidation.  At  a  higher  temperature, 
the  hydrogen  oxidation  rate  increases  [20,21].  Therefore,  op¬ 
timization  of  the  gas  mixture  temperature  at  the  PROX  reactor 
inlet  is  essential  to  realize  the  enhanced  catalyst  selectivity  to 
CO  oxidation.  The  inlet  oxygen-to-carbon  monoxide  molar 
ratio  is  an  important  parameter  in  the  CO  oxidation  process  in 
the  preferential  oxidation  reactor.  Fig.  13  shows  the  effect  of 
the  inlet  oxygen-to-carbon  monoxide  molar  ratio  on  the  ac¬ 
tual  CO  conversion  versus  the  catalyst  weight  profile.  From 
the  plots,  it  is  apparent  that  a  lesser  amount  of  the  catalyst  is 
needed  at  a  higher  oxygen-to-carbon  monoxide  molar  ratio 
than  that  at  a  lower  oxygen-to-carbon  monoxide  molar  ratio 
to  achieve  the  same  actual  CO  fractional  conversion.  This 
suggests  that  the  PROX  reactor  should  be  operated  at  a  value 
higher  than  the  stoichiometric  value  of  0.5  of  the  oxygen  to 
carbon  monoxide  molar  ratio  to  reduce  the  amount  of  catalyst 
needed.  However,  the  selected  oxygen  to  carbon  monoxide 
molar  ratio  should  be  such  that  the  hydrogen  (desired  fuel-cell 
fuel)  oxidation  is  still  insignificant  relative  to  the  CO  oxida- 


Fig.  14.  Effect  of  total  inlet  molar  flow  rate  on  CO  conversion  in  the  PROX 
reactor. 
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tion  process  on  the  Pt-on-Al203  catalyst  in  the  PROX  reactor. 
Fig.  14  shows  the  effect  of  the  total  reactant  molar  feed  rate 
to  the  reformer  on  the  actual  CO  conversion  versus  catalyst 
weight  profile  for  the  four  flow  rates  given  in  Table  2.  It  is 
quite  obvious  that  a  larger  amount  of  catalyst  is  needed  for  a 
larger  total  reactant  molar  feed  rate  to  the  reformer  for  a  fixed 
actual  CO  fractional  conversion.  This  is  so  because  the  PROX 
reactor  is  then  required  to  oxidize  a  larger  amount  of  CO. 

5.  Conclusions 

(a)  A  preliminary  mathematical  model,  based  on  the  prin¬ 
ciples  of  thermodynamics  and  chemical  reaction  engi¬ 
neering,  has  been  developed  to  predict  the  performance 
of  a  reformer  for  the  production  of  hydrogen-rich  gas 
mixture  from  a  hydrocarbon  fuel,  such  as  ft-decane  (here 
used  as  a  model  compound  for  the  JP-10  (CioFl^)  fuel), 
for  the  fuel  cell  applications.  Although  the  model  was 
originally  developed  for  the  carbon-to-oxygen  atom  ra¬ 
tio  of  one  with  ^-decane  as  the  fuel  and  air  as  the  oxidant, 
the  model  can  be  adapted  to  the  non- stoichiometric  feeds 
of  a  fuel  and  air  for  the  partial  oxidation  of  the  fuel  in 
the  Rh-on-y-Al203  catalytic  partial  oxidation  reactor  to 
produce  syngas  for  its  further  treatment  in  the  water-gas 
shift  and  preferential  oxidation  reactors. 

(b)  A  computer  program,  based  on  the  performance  model 
equations,  was  developed  and  implemented  to  numer¬ 
ically  simulate  the  performance  of  the  fuel  processor 
shown  in  Fig.  1. 

(c)  The  relation  between  the  predicted  exit  gas  mixture  tem¬ 
perature  and  the  inlet  temperature  of  the  gas  mixture  en¬ 
tering  the  CPO  reactor  may  be  approximated  by  a  linear 
relation. 

(d)  The  computed  performance  data  for  the  HTWGS  and 
LTWGS  reactors  show  similar  trends.  The  gas  mixture 
temperature  and  the  actual  CO  conversion  increase,  with 
an  increase  in  the  catalyst  weight,  to  their  equilibrium 
state  values.  Optimization  of  the  gas  mixture  temperature 
at  the  inlet  of  the  HTWGSR  as  well  as  of  the  LTWGSR 
is  important  to  avoid  sintering  of  the  catalyst  and  reduce 
its  weight  in  each  reactor  for  a  desired  overall  actual  CO 
conversion  in  these  reactors.  This  also  has  the  effect  on 
the  amount  of  water  needed  for  injection  into  the  wa¬ 
ter  injector/gas  cooler  unit.  One  should  make  it  certain 
that  the  amount  of  water  injected  is  sufficient  to  convert 
the  entire  amount  of  CO  in  the  gas  mixture  entering  the 
HTWGSR-external  cooler-LTWGSR  module. 

(e)  The  reaction  rate  for  the  CO  oxidation  in  the  PROX  reac¬ 
tor,  as  predicted  by  the  empirical  equation  (19),  shows  an 
interesting  behavior.  It  changes  from  a  low  to  a  relatively 
high  value  as  the  actual  CO  conversion  approaches  unity. 
One  may  explain  this  in  terms  of  the  competition  be¬ 
tween  carbon  monoxide,  hydrogen,  and  oxygen  for  their 
chemisorption  at  the  catalytic  active  sites  for  the  overall 
CO  oxidation  process. 


6.  Suggestions 

(a)  In  the  performance  model  presented  in  this  paper,  the 
assumption  of  complete  conversion  of  a  fuel,  such  as  n- 
decane,  and  oxidant  (air)  at  the  carbon-to-oxygen  atom 
ratio  of  1  was  used  due  to  the  non-availability  of  the 
chemical  kinetics  rate  equation  for  the  partial  oxidation 
reaction  on  Rh-on-Al203  catalyst  in  the  catalytic  partial 
oxidation  (CPO)  reactor.  It  is  suggested  that  first  the  ki¬ 
netic  rate  law  equation  for  the  partial  oxidation  of  a  fuel, 
for  example,  ft-decane  (ft-Ciotfe)  or  JP-10  (CioH^)  be 
developed  through  experimental  work.  Then,  this  rate 
law  equation  should  be  incorporated  into  the  current  re¬ 
former  performance  model  to  predict  the  CPO  reactor 
performance  behavior  as  well  as  that  of  the  overall  re¬ 
former  system  for  different  carbon-to-oxygen  atom  ratios 
or  for  different  fuel  and  air  feed  rates. 

(b)  It  is  suggested  that  a  small-scale  fuel  reformer,  based  on 
the  scheme  presented  in  Fig.  1,  be  designed  and  built  so 
that  the  experimental  program  can  be  initiated  to  acquire 
the  experimental  data  for  the  improvement/validation  of 
the  reformer  performance  model  presented  here. 
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